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Abstract. - We consider tunneling in a hybrid system consisting of a superconductor with 
two or more probe electrodes which can be either normal metals or polarized ferromagnets. In 
particular we study transport at subgap voltages and temperatures. Besides Andreev pair tun- 
neling at each contact, in multi-probe structures subgap transport involves additional channels, 
which are due to coherent propagation of two particles (electrons or holes), each originating 
from a different probe electrode. The relevant processes are electron cotunneling through the 
superconductor and conversion into a Cooper pair of two electrons stemming from different 
probes. These processes are non-local and decay when the distance between the pair of in- 
volved contacts is larger than the superconducting coherence length. The conductance matrix 
of a three terminal hybrid structure is calculated. The multi-probe processes enhance the con- 
ductance of each contact. If the contacts are magnetically polarized the contribution of the 
various conduction channels can be separately detected. 



Superconductor-normal metal (SN) contacts at mesoscopic scale are of primary importance 
in view of the interplay between coherence effects in the metal and intrinsic coherence of the 
superconducting condensate which is probed by Andreev reflection [|l],||]. The situation be- 
comes even more interesting with ferromagnetic metals (F). The subgap Andreev conductance 
tends to be hindered by magnetic polarization |^|-^ of the F electrode, since not every electron 
from the spin up band can find a spin down partner to be converted in a Cooper pair. Recent 
experiments [^,0] suggesting a proximity effect in SF junctions are still under debate. More 
generally, multi-terminal NS structures offer the possibility of manipulating phase coherent 
transport . An example is the Andreev interferometer consisting of a mesoscopic N sample, 
made of two arms connected to two S "mirrors" on one side and to a single reservoir on the 
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other, where electron propagation m N is sensitive to the difference of the phase of the two 
superconductors. In general the two SN contacts are separated by a distance larger than the 
superconducting coherence length ^ = ?iUi?/7rA, A being the superconducting energy gap. 
A dual configuration, consisting in a superconductor connected to two N probe electrodes 
separated by a distance smaller than ^ (see Fig l.a), with two independent reservoirs, was 
considered in Ref. It has been proposed that correlations between the N probes could 

be established across the superconductor, by a process where two electrons each originating 
from a different N electrode, are converted in a Cooper pair. Recently, it has been shown that 
such correlations could be used to build entangled states of electrons If the probe 

electrodes are ferromagnetic, the conductance due to these Crossed Andreev (CA) processes 
is sensitive to the relative magnetic polarization, being maximal for opposite polarization of 
the F probes 0. 

In the present Letter, we investigate in detail the dependence of the conductance on the 
distance between the probes and on their spin polarization. We consider a three-terminal 
device A/S/B (see Fig. la) where S is an s-wave superconductor and the electrodes A and 
B can be either normal metals or ferromagnets. We study the linear conductance matrix at 
subgap temperature and voltages (the subgap regime T, eVi ^ A is assumed in what follows). 
Then no single electron channel is left and all the relevant channels involve simultaneous (on 
a time ^ ?i/A) tunneling of two electrons or holes (Fig. lb) : (a) single-contact Andreev 
reflection {2 A and 2B), where two electrons, both originating from the same electrode, A or 
B, are converted in a Cooper pair giving rise to the currents G2aVa and G2bVb', (b) CA 
processes, where two electrons, each originating from a different A/ B electrode are converted 
into a Cooper pair, with associated current Gca{Va + Vb) from each A/B probe to S; (c) 



cotunneling (EC), which is easily visualized in the tunneling limit |13| as processes in which, 
for instance, an electron from A tunnels to B via a virtual state in S; EC processes yield a 
current from A to B, Gec{Vb — Va)- Then the subgap currents can be presented in a matrix 
form 



\ ( G2A + Gca + Gec Gca — Gec \ i 
Ib I \ Gca — Gec G2B + Gca + Gec / \ Vb 



(1) 



The multi-contact processes, CA and EC, increase the "diagonal" conductances diA/dVA, diB/dVB 
and give rise to off-diagonal terms dlA/dVB,dlB/dVA, i.e. the current at probe A (B) de- 
pends also on the voltage at probe B (A) . The generalization of Eq. (|l|) to more complicated 
structures (see Fig.l.c) is straightforward. 

For illustrative purposes we calculate the conductance matrix for contacts being tunnel 
junctions. The single-junction conductances G2A and G2B were calculated in Ref. ||l^. Here 
we find that EC and CA conductances depend on the relative position R of the contacts, 
vanishing exponentially for \R\ 3> ^. For multichannel tunnel junctions they are found to 
depend on propagation in 5* in a way which is sensitive to the geometry of the sample. This is 
due to interference effects between different channels, which also play a role in determining the 
single-junction NS conductances [Q. In multichannel junctions with normal A/B electrodes 
we find Gca = Gec, thus the off-diagonal terms in the conductance matrix Eq.(^ vanish. 
This symmetry is broken if A/B are polarized ferromagnets. Indeed Gec is suppressed if 
A and B have opposite polarization, because EC processes preserve the spin of the involved 
electron (in absence of any magnetic scattering whatsoever) and take advantage from parallel 
A/B polarization. On the other hand CA tunneling is suppressed for parallel A/B polarization 
since it is difficult to find two partner electrons with opposite spin to pair up. Then if 
magnetically polarized probes are used, off-diagonal conductances in Eq.(|l|) are non vanishing, 
due to the presence of unbalanced CA and EC processes, for inter-contact distance |^| <£^. 
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Fig. 1 - (a) Schematics of the three-terminal A/S/B device; (b) diagrammatic representation of the 
processes leading to subgap conductance: single-contact two-particle tunneling (2 A and 2B), elastic 
cotunneling (EC), which probes the normal Green's function (full lines) of S and "crossed Andreev" 
(CA), which probes the anomalous propagator (dashed line); (c) a simple experimental design for 
measuring the excess EC and CA currents: various electrodes {A,B,C',D,E) allow to probe different 
distances Ra. 



Cotunneling and Crossed Andreev tunneling rates. - We describe contacts in the tunnel- 
ing regime by standard tunneling Haniiltonians 

Hta = E t4 cldp, + T^f* d^cka ; Htb 4.c«- + ^if ^* (2) 

kpa- pqcr 

where and are matrix elements between single electron states fc G A, p G S' and q ^ B. 
Quasiparticle states in S are defined by the operators — Upa-d^^ — Vpad^p^„. 

We now classify processes which appear in perturbation theory in TYti- Single electron 
processes are absent in the subgap regime. In lowest nonvanishing order we have to con- 
sider only "elastic" processes, where the quasiparticle created in the intermediate state is 
destroyed "Inelastic" processes leave an excitation of energy > A in the superconduc- 

tor, so under subgap conditions they can be neglected. The relevant processes, represented 
in Fig. lb, are single-contact Andreev reflection in the tunneling limit {2 A and 2B) [|l^, CA 
tunneling, and "elastic" cotunneling [ p^ , where the same electron tunnels at once from A to 
B {B to A). One notices that EC and CA are non-local probes of electron propagation in 
the superconductor ^ . This is apparent if one considers the gedanken case of single channel 
tunnel junctions, where the size of the contacts is d^, ds ~ Ai?, the Fermi length. EC probes 
the "normal" propagator in the superconductor whereas CA processes probe the anomalous 
propagator (see Fig. lb) and both processes are exponentially suppressed for 

We now turn to the calculation of the spin dependent tunneling rates, using Fermi's golden 
rule. We consider electrodes A and B which can be magnetically polarized (parallel or an- 
tiparallel) along the same axis, resulting in different spectral properties. We neglect here the 
influence of ferromagnetic electrodes on the superconductor's spectrum. This a reasonable 
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assumption for small tunneling contacts on a massive superconductor. In this case, the prox- 
imity effect should not be too serious, e. g. superconductivity is not destroyed in the contact 
vicinity. This is in contrast with Ref. JTst where a sandwich geometry is considered, and where 
conduction above the gap is studied. 

By proceeding along the lines of Ref. we write the EC rate Aa Ba as 



r 



A- 



27r 



jdede'dQdC 5{e - e') f{e - cVa) [1 - f{e' ~ cVb)] 

F£c(C,e)F£c(C',e') S|;c(e-eV:i,e'-el/s,C,C') (3) 



where /(e) is the Fermi function, F£;c(^,e) = ((■ + + — contains information on 
energies of virtual quasiparticles in S. The rate for CA processes, {Aa, B — cr) — *■ 5, is given 
by a similar expression 

r^B^s = ^ Jdede'dCdC S{e + e') f(e-eVA)f{e'-eVB) 

FcAiC,£)FcA{C,e') E^cAi£-eVA,£'-eVB,CX') (4) 

where Fcyi(C, e) = A/(C^-|-A^ — e^). Information about propagation in the specific geometry 
are contained in the functions e', C')- Here we give the explicit expression for planar 
uniform tunnel junctions and local tunneling, T(r, r') — T 6(z) 5{f—r'). We moreover consider 
ballistic propagation of plane wave states (generalization to diffusive conductors can be made 
following Ref. 10). In this case the functions 'E.ec and can be expressed as 

S|;c(e,e',C,C') = \TaTb\^ I dridr2 ! dr^dr^ Jl(12, e) Jg(31, C) Jg(24, C') Jg(43, e') (5) 

J A J B 

S£^(£,£',C,C') = \TaTb? fdndr2 fdrsdn JX(12, e) Kg(31, C) Kg"(24, C') V(43, e') 

J A J B 

where the spectral functions are defined as, e.g. JJ(12,w) = J^(ri,r2,w) = Tlk^i'^ ~ 
Eko-) V'k<T(ri)Vic.(r2), Kg(31,a;) = KX(r3,ri,u;) = ^('^~^ka) V'-k-<T(r3)V'k<T(ri). The space 
integrals in (0) run on the contact surfaces. The diagrammatic representation is given in 
Fig. l.b. 

At low temperature and voltages the main contribution to the rates Eqs.(^,P is due to 
electrons close to the Fermi level in the A/ B probes, e = ±£' « 0, and the leading dependence 
of the rates on the voltages comes from the Fermi functions. One can let S(£,£',C, C') ~ 
S(0,0, C, C') and perform the e and e' integrations in Eqs.(^,^. At T = this gives I'^q = 
eV^^B = G^Eci^B - Va) and I^a = SeF^^^^ = 2Gg,^(VA + Vb), which define the spin- 
dependent conductances. 

It is instructive to consider first single channel junctions {dA, dB ^ Af)- The conductances 
are calculated by putting ri — r2, = r4, ri — = R in Eq.(|^). By performing the C and 
integrations the result is obtained 

/ Gpc \ 2TT^e^ Q-^RK ( N?,(0) cos^ (ksR) \ 

\G-cA I ^ \ sy > AV ) (^^^)2 ^^.(0) s\n\ksR) J 

where ks and Ns{0) are the Fermi wavevector and the normal state density of states of 
the superconductor. Propagation in the superconductor is characterized by factors which 
depend on |^|, being periodic with period ir/ks and being suppressed for |^| ^ The 
magnetic polarization of A/B electrodes enters only via the spin-dependent density of states 
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7V^^(0), assuming that -/V^^(O) = N^^ g{±^h) where h is the exchange field, Nj^ g is the 
density of states in absence of magnetism and =F sign stands for spin a (anti) parallel to the 
magnetization. If the density of states for minority spin can be neglected, then for parallel 
(antiparallel) polarized contacts only EC (CA) processes will yield a conductance. 

For the more realistic case of multichannel tunnel junctions, interference between different 
channels has to be considered. This has been studied in Ref. for single contact processes 
(^1 and A2) in a NS tunnel junction, which probe propagation in the normal electrode (it 
is in general diffusive and depends on the geometry). As for EC and CA, we notice that 
in clean superconductors the single particle propagator is rapidly oscillating (~ ks), so it 
would average out on a length ^ <C 1/ks in multichannel junctions. This is not the case for 
two-particle propagators involved in EC and CA processes, as we show explicitly below by 
considering a specific geometry. 

At this stage one could assume that the overall conductance is given by independent 
single-channel contributions, start from Eq.(|^) and argue that the factors cos^ (ksR) and 
sm^{ksR) are averaged over distances ~ dA, ds- The resulting expressions for the EC and CA 
conductances would be nearly identical yielding, in the special case of normal A/B electrodes, 
Gec = GcA- This conclusion turns out to be correct, even if the actual expression for the 
conductances involves interference between different channels. To show that we start the 
calculation from Eq.(j^), accounting also for different spin-dependent Fermi wavevectors fc^ 
(fc^) of the A (B) electrode and kg. We perform the coordinate integrations, and at this 
stage it becomes apparent that the terms leading to the special dependence cos^ {ksR) and 
sin^(fcsi?) drop out. Because of interference between different channels, the result depends 
on the geometry and on the mismatch between the spin-dependent Fermi wavevectors. The 
simplest case is a geometry with the two junctions belonging to the same plane (see Fig. la), 
where the results depend on the distance R > dA^ds 



G'hA J 8e2 \^ ^b'G^'^ / 



(7) 



Here the factors JFJ g contain information on the geometry and g are the spin-dependent 
one-electron conductances for each junction (S being in the normal state), for instance 



« Z:^ 7V-(0)iV5(0) H<,y^dA) (8) 

where ka,b = {k'A bI^sY^"^ and S is the area of the junction. Here the single junction 
geometry factor is given by J^{K?,y) = 2it ^^{dx/x) sin(Ka;) sin(x/K) and the two-junction 
factors in Eq.(0) are, e.g., Ta = [^(k^, Vfesfcl^A)]"^ Tie J^'^dO [1 - e»'=sdA(Ki+cos e)]/(^2^ ^ 
cos 6) for the geometry we consider. It is important to point out only some general property. 
For K = 1 the factor y) (x\ny depends weakly on the reduced size y; the factors J-a,b 
are even more weakly dependent on y and substantially of order one. A slight asymmetry 
K ^ 1 makes all the factors independent on the size of the junctions, if dA,dB are large 
enough, y\K — k~^\ ^ 1. Still Ta,b are of order one, so EC and CA processes determine an 
appreciable conductance. 

Discussion. - We can now discuss the full conductance matrix in equation (|l]) by defining 
the total EC and CA conductivities Gec = Gl;^ + G^^ and Gca = Gg,^ + G^^. The EC 
and CA conductances appear both in the diagonal and in the off-diagonal conductance matrix 
elements in Eq.dll). Let us first consider the case of non magnetic contacts, where we can drop 
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the spin dependence. For multichannel contacts Eq.(^ shows that Gec — Gca , so the off 
diagonal conductances vanish. Coherent tunneling processes involving two distant contacts 
enter only the diagonal terms, and provide as an extra contribution with respect to the 
standard Andreev conductances G2A and G2B. The extra current depends on the distance R 
between the two contacts. A simple setup where the R dependence of the extra current can be 
studied is shown in Fig.lc (alternatively one may use a STM tip as a mobile contact). Another 
signature of CA and EC processes can be found if one considers contacts of very different 
transparency, say \Ta\ ^ \Tb\- In this case G2A oc \Ta\^ dominates Gec+Gca oc iT^plTep, 
which itself is much larger than G2B ~ I^^bI*- Thus the conductance at the less transparent 
probe (lower right diagonal element in Eq.(j^)) is given by Gec + Gqa, so it is essentially due 
to two-contact processes: if we bias contact B a current will flow because of correlations of 
superconductive nature with contact A. li R < f, the crossed conductances are still affected 
by the factor (ksR)^^ which can be very small ||l2| . This problem can be partially overcome 
for instance by choosing one contact A to be of size >> ^ or even a semi-infinite interface. 

Let us now consider the case of spin polarized probes. As it is apparent from Eq.(^ 
Gec 7^ Gca so the off diagonal elements in the conductance matrix, Eq. (^ , are finite and the 
current in one contact can be manipulated by the voltage bias of the other contact. The sign 
of this effect depends on the mutual polarization of the electrodes. This generalizes the result 
of Ref. [|l0|. More in detail, spin polarization enters in two ways in the result (|^,0): first, in 
the spin-dependent densities of states; second, in the shift of the Fermi momenta fc^ 3, which 
modifies the factors in Eqs.(0). To fix the ideas, let us for simplicity neglect the latter and 
concentrate on the effect of the density of states. Defining the contact polarizations Pa,b = 

Tpr§-, one has simply that Gec is proportional to (1 + PaPb) and Gca to (1 — PaPb)- 

Therefore the off-diagonal conductance is roughly proportional to {—PaPb)- This shows a 
striking consequence of the competition between cotunneling and crossed Andreev processes, 
via their spin-dependence : non only the amplitude, but also the sign of the conductances 
can be controlled by spin polarizations. In the extreme case of parallel complete polarizations 
the only possible process is cotunneling, with I a = —Ib, while for antiparallel polarization 
crossed Andreev tunneling prevails, with I a = Ib- 

So far we have discussed the zero-temperature case. Direct generalization to finite tem- 
perature leads formally to a divergence of the tunneling rate. It is due to the finite, though 
very small (oc e~^/^), probability of exciting an electron from A{B) to a quasiparticle state 
in S, and to the divergence of the quasiparticle density of states in S. The divergence in the 
rates disappears if the latter is rounded off at A. The EC and CA conductances acquire 
an additional contribution (x e~^/^ ln(A/r) where F is a scale related to the mechanism of 
broadening of the quasiparticle levels in the superconductor. 

In the present Letter we have demonstrated the non-local character of cotunneling and 
Andreev reflections on a superconductor and we have also studied the role of magnetic po- 
larization. We have discussed possible schemes to detect these effects in devices with three 
or more terminal. Devices with high transparency contacts [ p^ are also promising for exper- 
iments. Further theoretical analysis would require the self-consistent analysis of the mutual 
effects of superconductivity, diffusive propagation and ferromagnetism in the hybrid system. 
For both high and low transparency contacts propagation in the specific geometry has to be 
taken into account. 



* * * 
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